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Abstract 



We prove the existence of reaction-diffusion traveling fronts in 
mean zero space-time periodic shear flows for nonnegative reactions 
including the classical KPP (Kolmogorov-Petrovsky-Piskunov) non- 
I> | linearity. For the KPP nonlinearity, the minimal front speed is char- 

acterized by a variational principle involving the principal eigenvalue 
(f-} . of a space-time periodic parabolic operator. Analysis of the variational 

r***'" \ principle shows that adding a mean-zero space time periodic shear flow 

to an existing mean zero space-periodic shear flow leads to speed en- 
. hancement. Computation of KPP minimal speeds is performed based 

on the variational principle and a spectrally accurate discretization of 
the principal eigenvalue problem. It shows that the enhancement is 
monotone decreasing in temporal shear frequency, and that the total 
enhancement from pure reaction-diffusion obeys quadratic and linear 
laws at small and large shear amplitudes. 
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1 Introduction 



Front propagation in heterogeneous fluid flows is a fundamental issue in com- 
bustion science, and has been an active research area for decades (see jl] , [H] , 
[Oj, j20ji |22], j2Hj [2Z1 and references therein). A fascinating phenomenon is 
that the large time front speed depends strongly on the flow structures. Front 
dyanmics and speed enhancement have been studied mathematically for var- 
ious flow patterns by analysis of the proto-type models, i.e. the reaction- 
diffusion-advection equations (see 00011011111113113111111111221123123 
and references therein) . In ^3] , ^7] , fronts in space-time periodic shear flows 
(defined below) have been studied both analytically and numerically. This 
paper is a continuation of |T7j to establish front existence and variational 
principle on front speeds. The variational principle will be utilized both for 
analysis and computation to gain qualitative and quantitative understanding 
of front speeds. 

Let us consider the reaction-diflusion-advection equation: 

u t = A x u + B- Va.it + /(«), (1.1) 

where u = u(x,t), x G R n , n > 2, t G R; A x and V x are the standard 
Lapliacian and gradient operators in x. To specify the vector field B in R n , 
let us write x = (xi,y), with y G R n ~ l . Then B is a smooth, time-dependent 
shear flow B = (b(y, t), 0, . . . , 0) where b(y, t) is a scalar function with period 
equal to one in both y and t, and mean equal to zero over the period cell 
(denoted by T n ): 

I b{y,t)dydt = 0. (1.2) 
To exclude degeneracy, we assume that: 

/ \V y b(y,t)\ 2 dydt^0. (1.3) 

We shall consider nonlinearity / : [0, 1] — > R to be smooth and of KPP 
type: 

/(0) = f{l) = 0, f(u) > 0, V u e (0, 1), /'(l) < 0. (1.4) 

The definition (jl.4j) includes functions such as f(u) = u m (l — u), m > 1, 
integer; and the Arrhenius-type nonlinearity f(u) = e~ E l u {\ — u), E > 0, 
arising in combustion [23- The classical (strict) KPP nonlinearity [23 obeys 
the additional inequality f{u) < uf'(0), Vw G (0, 1), e.g. f{u) = u{\ — u). 
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We are interested in traveling fronts moving in the direction along x\ and 
of the form: 



u = U [x\ — ct, y, t) = U(s, y,r), s = x\ — ct G R, r = t, (1.5) 
with U being 1-periodic in both y and r and satisfying conditions at s = ±00: 

lim U(s, y, r) — 0, lim U(s, y, r) — 1, (1.6) 

s— >— 00 s— >oo 

uniformly in ?/ and r. After substitution, the problem becomes: 

A a , y U + (c + b(y, r)) U 9 - U T + f(U) = 0, (1.7) 
lim U(s, y, r) = 0, lim U(s, y, r) = 1, U has period 1 in (y, r), 

s— >— 00 s— »oo 

which has a parabolic equation on RxT n with both the function {7 and the 
constant c unknown, or a nonlinear eigenvalue problem. 

The approach to front existence is to start with the front solutions to (jl-H) 
with f(u) being a so called combustion nonlinearity with cut-off: fe(u) = 
for u G [0,9], 6 G (0, 1), f e (u) > for u G (0, 1). = and //(l) < 0. 

In ^7], the present authors outlined the existence proof of such combustion 
fronts using the method of continuation |24. Here we provide the details, 
with the front speed estimates drawn from The existence of KPP type 
fronts traveling at minimal speeds follows as suitable limits of the combustion 
fronts as 9 — > 0. Such a construction method of KPP type fronts is due to 
[3] on fronts in spatially dependent shears. 

In case of strict KPP nonlinearity, we show that the minimal speed de- 
noted by c* is given by the variational principle: 

c* = - inf (1.8) 
A>0 A 

where \i = //(A) is the principal eigenvalue of the periodic parabolic operator 

L x = A y - -d T • +[A 2 + A b(y, t) + f (0)], (y, r) G T n . 

The variational principle (jl.8|) allows us to deduce enhancement proper- 
ties of c* from knowledge of /1. Moreover, (|1.8J) provides the most efficient 
way of computing c* without solving the original time dependent PDE (jl.lj) . 
We shall first compute fi by a spectrally accurate method, then minimize. 
With more accuracy in much shorter time, the minimal KPP front speeds 
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are computed. The quadratic and linear speed enhancement laws from the 
pure reaction-diffusion speed (without advection), as well as monotonicity of 
enhancement in shear temporal frequency are demonstrated, in agreement 
with [T7j based on other methods. In particular, both analysis and computa- 
tion from the variational principle show that adding a mean zero space-time 
shear to a spatial mean zero shear alway generates additional enhancement. 

The paper is organized as follows. In section 2, we present properties 
and existence proof of combustion fronts. In section 3, we prove existence 
of KPP type fronts with minimal speeds. In section 4, we prove the vari- 
ational principle ()1.8|) and its qualitative properties. In section 5, we show 
computational results of KPP minimal speeds based on (|1.8|) . 

2 Combustion Fronts and Properties 

In this section, / is the combustion type nonlinearity with smooth cut-off. 
We prove: 

Theorem 2.1 (Existence and Uniqueness) Let f be of combustion type. 
Then there exists a classical solution (U, c) to the problem \1. 7| ). Moreover, 
this solution {U,c) is unique up to translation in s. That is, if (Ui,ci) and 
(JJiiCz) Qjre both solutions, then c\ = C2 and Ui(s,y,r) = L^s + X,y,r) for 
some A G R. 

The proof will be preceded by a list of properties of front solutions. 

2.1 Strong Maximum Principle and Comparison 

As in jS] and [21] , one of the most important tools in our analysis is the strong 
maximum principle. We denote by L c the linear part of the operator in ([1.7)1 . 
This operator is parabolic, so the usual parabolic maximum principle and its 
corollaries hold. Specifically, if U attains its minimum at some finite point 
(sq, yo, r ), then U is constant for r < t . The functions on which L c operates 
are periodic in the r variable (time), so the constant extends to all r G R. 
This gives us a strong maximum principle for L c that is analogous to the 
maximum principle for elliptic operators: 

Proposition 2.1 (Strong Maximum Principle) Suppose U G C 2 and 
L c u = A StV U+ (c + b)U s - U T < 0. (2.1) 
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Suppose also that U attains its minimum at some finite point (sq, jjq, t ). 
ThenU is constant for all (s,y,r) G RxT n . IfH(s,y,r) < and L c u+Hu < 
for all (s, y, r) and u attains a nonpositive minimum at some finite point, 
then u is constant in Rx T n . Furthermore, if L c u + Hu < and u > in 
R x T n , then ifu(so, yo, To) = at some finite point, then u must be constant 
in R x T n , even if H < does not hold. 

Since solutions to (IT7I) satisfy LJJ < -f(U) < 0, it follows that < U < 
1. By the method of sliding domain [3], [21], the strong maximum principle 
implies the following two Theorems. 

Theorem 2.2 (Monotonicity) Let the nonlinearity f be combustion type. 
Then the solution U of \1.T\) is strictly monotone in s: U s > for all 
(s,y,r)eRxT n . 

Theorem 2.3 (Comparison Theorem) Let f E C 1 ([0,1]) be a combus- 
tion type nonlinearity. Let /3(y,r) and /3(y,r) be smooth and periodic with 
(3 > ft for all (y, r) e T n . Assume that u and u G ^R 1 x T n ) are two 
functions that satisfy 

Lu + f{u) = A. %y u + ft(y, t)u s - u T + f(u) > 0, 

Lu + f(u) = A S:V u + ft(y,r)u s -u T + f(u)<0, (2.2) 

and the conditions at s infinities: 

lim u = and lim u — 1 

s— oo s^+oo 

lim u = and lim it = 1 

s— >— oo s^+oo 

Additionally, assume that at least one of the following holds: 

ft = ft OR u s > OR u s >0. (2.3) 

T/ien there exists X e R such u(s,y,r) = u(s + X,y,r) for all (s,y,r). 
Furthermore, there exists a point (yo, r ) G T n such that ft(yo, To) = ft(yo, To), 
so i/iai ft > ft cannot hold for all (y, r) G T n . 

Proof: The proof of Theorem 2.3 uses the sliding domain method. Define 
the function ui\ = u(s + A, y, r) — u(s, y, r). The key is that the inequalities 
()2.2j) are invariant under translation in the s direction. For all A G R, 

Lw x + H x w x <{ft- ft)u s (s + A) (2.4) 



5 



where 

H x = H x (s, y,r)= [ f(&(s + A) + (1 - £)u) d£. (2.5) 
Jo 

Assuming either ft = ft or u s > is satisfied, the right hand side of (|2.4|) is 
non-positive. In the event that u s > 0, we consider instead 

Lw x + H x w x <(ft- ft)u s (s + A), (2.6) 

so that again the right side is non-positive. In either case, w x satisfies 

lim w x = and lim w x = 0. 

Using the strong maximum principle for L (or L), we can show that for A 
sufficiently large, w x > everywhere in R x T n . Then, for \i defined by 

fi = inf {A \ w x > in Rx T n , V A > A } , 

it can be shown that = 0. In other words, we can translate u to the 
left so that u > u. Then translating back to the right until they touch, we 
have u > u with equality holding at a finite point, and the strong maximum 
principle implies the result. If ft > (3, the right hand side of inequality ()2.4|) 
becomes strictly less than zero for all A. Plugging in w^, we then obtain a 
contradiction, implying ft(y ,T ) = (3(yo,To) at some point. The proof relies 
on the fact that f(u) vanishes for u < 9 and that /'(l) < 0, so that H x < 
in the region where \s\ is sufficiently large. Note that this is not necessarily 
true for the KPP nonlinearity because in that case f(u) > for all u G (0, 1). 

Theorem 2.2 follows from Theorem 2.3. We take u = u and ft = ft and 
show that /i = 0, hence, u s > 0. Differentiating equation (jl.7|) with respect 
to s, we obtain an equation for u s . The strong maximum principle and the 
limits as s — > ±oo then imply strict inequality: u s > 0. 

2.2 Exponential Decay 

As with the problems considered in [21] and j3] , we can obtain precise infor- 
mation about the decay rate of solutions through related eigenvalue problems. 
Specifically, solutions of (jl.7|) must decay exponentially fast to the limits 
and 1 as s — > ±oo, as described in the following theorem: 



6 



Theorem 2.4 (Exponential Decay) If (U,c) is a classical solutions of 
\1. 71 ) with combustion nonlinearity, then there exist constants, C\ > 0, C 2 > 
; Ai > ; A2 < 0, and positive functions $i(y,r), $ 2 (y,r) G C 2 {T n ) such 
that 

U(s,y,r) < Cie Al '$i(y,r), s < Sl 
l-C/(s,y,r) < C 2 e A2S $ 2 (y,r), s > s 2 (2.7) 

Proof: This result is analogous to Proposition 1.1 in [21], the only significant 
difference being that the related eigenvalue problem is parabolic. If Si < 
is such that U(s, y,r) < 6 for s < Sj, then we have 

L C C/ = for s < si. (2.8) 

To compare U with an exponentially decaying function, we want to solve 
L c w = for some function w of the form w = Cie^ s &(y, r) and $ > 0. Such 
a function must satisfy 

e- Xs L c w = - $ r + [A 2 + A(c + b(y, r))]$ = 0. (2.9) 

We will call p(A) G R the principal eigenvalue of this operator if for some 
$ > 0, 

A y $ - $ T + [A 2 + A(c + 6)]$ = p(A)$. (2.10) 
Thus, we want to find some A > such that p(A) = 0. 

Lemma 2.1 For each A G R there exists a principal eigenvalue p(X) and a 
principal eigenfunction $ > solving A2.1U\) . 

Proof: This is a periodic-parabolic eigenvalue problem, and again we see the 
importance of the periodicity of the time dependence that leads to a strong 
maximum principle. Problems of this type have been treated in more detail 
in and [12]. Let 

E = C 2+a ' 1+a/2 (T n ) and F = C a ' a/2 (T n ). (2.11) 

Let M : E — > F be the operator defined by, 

M$ = $ T - Aj,$- [A 2 + A(c + &)]$. (2.12) 
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or M$ = (- — d)$>, so that p(X) = - — d and $ are the desired eigenvalue 



For fixed A we can pick d > sufficiently large so that the operator (M + d)<& 
is invertible (see Section 2 of [12 J. By parabolic Schauder estimates, the 
map T = %{M + d)~ l : F — > F is compact, where i is the embedding of 
E into F. Since the strong maximum principle implies positivity of the 
operator T, the Krein-Rutman theorem then implies the existence of a simple 
eigenvalue \x with positive eigenfunction <3> for operator T. Hence T$ = p<3>, 

i - so tW ~ m - 1 
and eigenfunction for M. 

Because of the periodicity in y and r, integrating the equation over T n 
gives us 

A 2 / <5>dydT + \i (c + b)<5> dy dr = p{\) <5>dydr, (2.13) 
Jt" Jt™ Jt™ 

and since $ > we can divide by f Tn $ to obtain 

P (A) = A 2 + Ac + Aip^. (2.14) 

The quotient on the right is bounded in absolute value by A||&||oo, so we see 
that as A — > oo, p(A) — > +oo. Setting A = it is also clear that we must 
have p(0) = with $|a=o= 1 (i-e. positive constant). Taking the derivative 
of equation (|2.1()j) with respect to A and setting A = and $|a=o= 1 we 
obtain 

A w * A -(* A ) T + (c + 6)=p A (0). (2.15) 
Averaging over y and r we get 

p A (0) = c, (2.16) 

since b(y, r) is mean zero and $ is periodic. This implies that 

p(A) = cA + 0(A 2 ). (2.17) 

If we integrate equation (|1.7|) over all s, y, r, then we see because of period- 
icity, integration by parts gives us 



c=- f(U)<0. (2.18) 

JRxT n 

Thus, for A small and positive, we must have p(A) < 0. Since p(A) — ► +oo 
as A — > oo, there must be some Ai > such that p(\\) = 0. The principal 
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eigenfunction, $ 1( corresponding to eigenvalue p(\i) is strictly positive, and 
w = e^ lS $i solves Lw = 0. 

To prove exponential convergence as s ^ +00, we perform the same 
analysis on the function V = 1 — U in the region [S2, +00) x T n , where S2 > 
is such that 

/'([/)< < 0, Vs>s 2 . (2.19) 

2.3 Existence of Traveling Front Solutions 

Now that we have have established the strong maximum principle for L c , the 
exponential behaviour of solutions as s — > ±00, and the existence of eigen- 
values for the related eigenvalue problem (|2.1(Jj) . the proof of Theorem 2.1 
follows from the reasoning in [21] with modification to account for the time- 
dependent component of operator. We work with the equation in the whole 
space RxT n and prove existence of solutions via the method of continuation. 

First, supposing that for some b = bo we have a solution (Uq, cq) to 
the problem 1)1.7]) . we find that there exist solutions corresponding to small 
perturbations: 

b s = b + 8b[ . 

We look for solutions in the same functions spaces as in [21] with the addition 
of the r dimension to the domain of the functions. As in [21] , we consider the 
linearized operator L Co v + f'(Uo)v and find that, because of monotonicity of 
solutions to 1)1.7]) . this linearized operator is Fredholm with one-dimensional 
kernel. Then, the existence of a solution (U s , c s ) corresponding to shear b s 
follows from a contraction mapping principle for S small. 

Next, to complete the continuation, we suppose there is some maximum 
value D > such that we can find solutions for all 5 G [0, D), and we want 
to obtain a solution for 5 = D, thus implying that no such maximum exists. 
To do so, we consider a sequence of solutions (U s , c 6 ) to the problems: 

A StV U 5 + (c s + b 5 )U 5 s - U 5 T + f(U 5 ) = 0, (2.20) 
lim U 5 = 1, lim U s = 0, U 5 has period 1 in (y, r). (2.21) 

We show that there is a subsequence (U S] ,c Sj ) converging to a solution in 
the case 5 = D. First, we need to bound the wave speeds c 6 in order to 
extract a convergent subsequence. To do this, we compare these speeds to 
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the wavespeed of the corresponding one-dimensional problem: 

v + kv + f(v) = 0, (2.22) 
lim v = 1, lim v = 0, 

s^oo s— >— oo 

with normalization 

v(0) = i (2.23) 

Lemma 2.2 Lei (f/ 5 ,^) sofoe iTOZ}) and let (v,k) solve <TQlj) . T/ien fc < 

C<5 + Halloo- 

Proof: Defining w(s, y, r) = t>(s) and w = U s , we have 

A S)2/ w + (3u s -u T + f(u) = 0, 
A s>y u + (3u s - u T + /(u) = 0, 

with /3(y, t) = k and /? = c 5 + & 5 (y, r). If k > c s + ||& 5 ||oo, then /3 > /3, so we 
can apply Theorem 2.3 to conclude that there must be some (yo, r ) G T n 
such that (3(yo, r) = /3(?/o, r o)- This implies k = c s + b 5 (yo, r ) < c 6 + ||& 5 ||oo- 
Now for 5 G [0,-D), Lemma 2.2 implies there exists constant < 
independent of 5 such that K < c s < 0. Hence, we can pick a subsequence 
{c Sn } that converges to some c* < 0. To have a solution ([/, c) with 5 = D, 
however, we must have c* < 0. The following lemma was proved in the 
appendix of [T7j : 

Lemma 2.3 (Bounding away from zero) If {(U s ,c s )} is a set of so- 
lutions to \2.20)) with 5 G [0,-D), then there is a constant 7 > independent 
of S such that 

c s < -7 < 0, V 5 G [0,-D). (2.24) 

Hence, c* < 0. 

Now, by standard interior parabolic estimates ([H], JH|), we can find a 
function U* and a subsequence such that on any compact set Q C 

Rx T n , U s i — > £7* uniformly with £/* solving 

A, iV ET + (c* + 6 D )[/ S * - t 7 * + /(f 7 *) = °- ( 2 - 25 ) 

We need to show that the limits (j2.21|) hold for U*. Since solutions 
are invariant with respect to translation in the s variable (that is, we can 
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translate them and they are still solutions), we can impose the normalization 
condition 

max U 5 (0,y,r) = 6 (2.26) 

on all solutions (U s ,c 5 ). By monotonicity in s, U 5 (s,y,r) < 6 for (s,y,r) G 
(—oo,0) x T n . Moreover, the limiting function U* must satisfy U* > for 
all (s,y,r). 

First, we show that lirn s _ > _ tX) U*(s, y, r) = 0. In the proof of Theorem 
2.4, we showed that for each 5 we have some C 5 , X s , <& s such that 

U s (s,y,r)<C s e^ s ^ s (y,r), s<0 

(here we take s\ to be because of the normalization condition). Because 
the eigenf unctions & s are unique up to multiplication by a constant, we may 
assume that L„ $ 5 = 1 for each <5. Since the wave speeds are bounded away 
from zero, c 5 < —7 < 0, the corresponding X s must be bounded away from 
zero. Supposing this were not the case, we take a sequence such that 

X s ' -> and $^ -> $ in C^T") with $ solving 

A y $ - $ r = 0. (2.27) 

The strong maximum principle implies that $0 = 1, which means that 
J Tn &<3?o = since b is mean zero. However, using (j2.14j) and j Tn $ Sj = 1, we 
have 

X Sj = -c s i - I > 7 - / (2.28) 

Letting j — > 00 we obtain a contradiction since the integral on the right 
hand side goes to zero. Thus, there exists a constant A* > such that 
X s > X* for all 5. Harnack's inequality implies that there is a constant 
M > independent of 5 such that 

inf $ 5 > M sup $ 5 > M. (2.29) 

Now chosing C* = 1/M, we have 

U 5 {s,y,r) < C*e Xs whenever s < 0, V (y,r). (2.30) 
This implies the desired result: lim s ^_oo U*(s, y, r) = 0. 
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Now we show that lim^oo U*(s, y, r) = 1. Since U* > 0, we can define 

lim U*(s,y,T) = ^(y,r), (2.31) 

s— >+oo 

and we want if) = 1. It is not hard to see that the function if) must solve the 
equation 

A y iP -ip r + f{if>) = (2.32) 

for (y,r) G T n . To see this, let h(y,r) be a smooth function on T n . Also, 
let £(s) : R — > [0, 1] be smooth with £ = for s < 1 and £ = 1 for s > 2. 
Multiply equation ()2.25|) by the test function 0(s, r) = h(y, r)£(s), average 
over [0, L] x T n : 



7 / / &s, v U*<]> + {c + b)V* a <j> + Ui<l> + f(U)4> = Q. (2.33) 



Parabolic estimates imply the boundedness of the s derivative of U*. As a 
result, integrating by parts and taking the limit as L — > 00 gives us the result 



i(;Ayh + if)h T + f{if))h = (2.34) 
for any test function h G C°°{T n ). Thus, if) solves 

A y ip - Vr + /(VO = (2.35) 

for (y, r) G T n . Now, the maximum principle implies that ■?/> must be con- 
stant. Hence, 

A y iP-iP T = = -f(iP), (2.36) 

so ip = 6 or ip = I, since /(VO > for if> G (0,1). If if) = 9, then by 
monotonicity, U* < 9 for all (s, y, r) so that U* satisfies 

Lu = A SjV U + (c + b D )U s -U T = -f(U) = 0. (2.37) 

Also, U* attains its maximum at some finite point (0, yo, Tq) due to the nor- 
malization condition. So, the strong maximum principle implies that 

U* = 9, (2.38) 

which contradicts the limiting behavior as s — > —00. Therefore, if) = 1 and 
we have shown that 

\im U*(s,y,r) = 1. (2.39) 

s— >oo 
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Now that we have a solution (U, c) to problem ([1.70 . the uniqueness of {7 
(up to translation in s) follows easily from the comparison theorem. Let (U, c) 
and (U, c) be two solutions to ([1.7(1 with c > c. Theorem 2.3 implies that we 
must have c = c and that for some constant A G R, U (s, y,r) = U(s + A, y, r) 
for all (s, y,r) eRx T n . 

Remark 2.1 Although periodicity in time plays a crucial role in the proof 
of these results, the periodicity of solutions in the y variable is not necessary. 
That is, if uo C i? n_1 is a bounded region with smooth boundary, our results 
generalize easily to treat traveling fronts in the cylinder (x, y) 6 {R x uj) when 
the shear is periodic in time with Neumann boundary conditions imposed in 
the y dimension: |^ = for y £ dcu. With the Hopf Lemma, the strong 
maximum principle of Proposition 2. 1 still holds, as well as the rest of the 
results of this section. The same is true of our results for the KPP equation 
described in following sections. 



3 Existence of KPP Type Fronts 

Using the results for the combustion case, we now construct a solution to the 
problem 

A Sty U + (c + b(y, r))U s -U T + f(U) = 0, (3.1) 
lim U = 1, lim U = 0, U has period 1 in (y, r), (3.2) 

s^oo s— >— oo 

where f(u) is KPP type nonlinearity ([1.4(1 . Since the shear b(y,r) will be 
the same throughout this section, we will use the notation L c to denote the 
operation 

L C U = A s>y U + (c + b(y, t))U. - U T . 

Theorem 3.1 (Existence) Let f 6 ^([0,1]) with f(u) > in (0,1), 
/(0) = f(l) = 0. Then there exists a classical solution (U*, c*) to the problem 

Proof: This proof extends the arguments in section 8 of [H] to the case with 
time-dependent coefficients. Following j3j, we let Xe( u ) £ [0, 1] be a smooth 
cutoff function satisfying 

X e(u) = if u < 9/2, 
Xe(u) = 1 if u > 9, 
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and xe' > Xe whenever 9' < 9. Each function fg(u) = f(u)xe{u) is a nonlin- 
earity of combustion type and fg — > / uniformly on [0, 1] as 9 — > 0. We will 
also assume that if 9' < 9, then for some point uq G (0,0), Xe'( u o) > Xo{ u o) 
and thus fg>(u ) > fg(u ). By Theorem 2.1, for each e (0, 1/2) there is a 
unique pair (u$,cg) solving 

L ce ue + fe(ug) = 0, (3.3) 
lim uq = 1, lim -u^ = 0, 

with normalization condition 

. 1 

maxu e {0,y,T) = -. 

We will show that the solutions (ug, eg) converge to a solution (u*, c*) of the 
KPP problem and that this speed c* is the minimum possible speed (minimal 
in absolute value) of any solution to ()3.1|) . First, we bound the sequence of 
wave speeds in order to extract a convergent subsequence. 

Lemma 3.1 The sequence of speeds eg is strictly decreasing as 9 — > 0. 
Moreover, the sequence is bounded, so that there exists c* < such that 
limg^o eg = c* . 

Proof: For the first claim, let 9i < 6\ < \ with corresponding solutions 
(wi,Ci) and (u 2 ,c 2 ): 

L Cl Ui + fg x {u x ) = 0, 

and 

L C2 U2 + fg 2 (u 2 ) = 0, 
L C2 u 2 + fg 1 (u 2 ) < 0. 

Arguing by contradiction, suppose c 2 > c\. Then for f3 = c 2 + b(y,r) and 
/3 = ci + b(y, r), we have ft > /3, so we can apply Theorem 2.3 to conclude 
that ft(yo, To) = ft(yo, To) for some point (y , r ) and that Ui(s, y, r) = u 2 (s + 
X,y,r) for some translation A G R. This implies that c\ = c 2 and that U\ 
must satisfy both 

L ci ui + f ei (ui) = 0, (3.4) 
L ci u x + fg 2 (ui) = 0. (3.5) 
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This cannot be true, however, since for some value of u%, /^(wi) < fg 2 (ui), 
by construction of the functions fg 1 ,fg 2 - Therefore c 2 < c\. 

As in [3], we obtain a lower bound on the sequence {eg} by a comparison 
with the wavespeeds of the corresponding 1-D problems: 

vg + k e v e + fe{ve) = 0, (3.6) 
lim vg = 1, lim vg = 0, 

s^oo s^—oo 

with normalization 

v,(0) = ~ (3.7) 

Using k = kg and / = fg, it follows from Lemma 2.2 that for each 9 e (0, 1/2), 
kg < eg + || & || oo - Thus, a lower bound on kg implies a lower bound on eg. 
Such a bound was derived in section 8 of 0. 

Since all the fg and ug are uniformly bounded, then as before, parabolic 
estimates imply the existence of a subsequence ug k converging uniformly on 
compact sets to a function U* that solves ()3.1j) . 

We must show, however, that the limits (j3.2j) are achieved. We know that 
for each 9, ug is monotone in s. This implies U* > 0. So we can define the 
limits 

^ + {y,T) = lim U*(s,y,r) and ip~ (y , t) = lim U*(s,y,r). (3.8) 

s— >oo s^— oo 

As in equation ()2.33|) . we multiply by a smooth test function and integrate 
by parts to show that that ip + and ip~ satisfy 

A^+-W + = -/(V + ), 

A^-C = -/(VO- (3-9) 

Integrating both equations over T n shows that f{ip + ) = and f{4>~) = 0. 
so that ip^ = const., i.e. or 1. Because of the normalization condition on 
the functions Ug, however, we must have 

max U*(0,y,r) = 1/2, (3.10) 

(y,r)eT" 

so that ip~ = and ip + = 1. Thus, the limits ()3.2|) as s — > ±oo are acheived, 
concluding the proof of Theorem 3.1. 

In the following sections, we will use c* to denote this unique speed ob- 
tained by the above construction of the solution (U*,c*). 
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Proposition 3.1 If (u,c) is any solutions to \3. 1\) and h3. fy) . then c < c* < 
0. 



Proof: Suppose (u,c) is a solution to ([3.1)1 . For any e > 0, we can pick 
6 > sufficiently small such that eg G (c*,c* + e), where (uq,c$) is the 
corresponding solution to the problem with combustion nonlinearity ()3.3|) . 
as in the construction of (u*, c*). Because fg < f, we have 



Arguing by contradiction, assume c > eg. We know that (ug) s > 0, so we 
can apply Theorem 2.3 with j3 = c + b(y, r) > eg + r) = (3 to conclude 
that c = eg. Hence c < eg < c* + e. Letting e — > 0, we conclude that c < c*. 

4 Variational Characterization of c* 

The existence results in the previous section hold for a general KPP nonlin- 
earity. If we impose the additional condition f(u) < uf'(0) for all u G [0, 1], 
we can obtain a variational characterization of the minimal wavespeed c* in 
terms of a related eigenvalue problem on T n , extending the characterizations 
given in [H] and j2] to this problem with a time-dependent shear. 

Let //(A) be the principal eigenvalue of the operator L x = A y ■ —5 T • + [A 2 + 
Xb(y, t) + f'(0)} with associated principal eigenfunction <p(y, r) > 0. That is, 



The variational characterization is given in the following theorem: 

Theorem 4.1 If the nonlinearity f is KPP-type and satisfies f{u) < uf'(0) 
for all u G [0, 1], then the minimal speed c* < is given by the variational 
formula 



L c u + f(u) = 0, 
L c u + f e {u) <0, 



and 



L C0 Ug + fgiug) = 0. 



L x <p = A y( p + [A 2 + Xb{y, t) + f'{0))<j> - <f> T = /x(A)0 



(4.1) 



c* = 




(4.2) 
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Proof: We define the value 



7* = - inf 4^- (4-3) 

A>0 A 



Clearly, 7* > —00. Let c < satisfy c < 7*. We show first that we must 
have c < c*, implying that 7* < c*. Since c < 7*, there exists Ao > such 
that 

-A c > /jl(Xq). 

Letting O = (f ) o(y, T ) > be the principal eigenfunction associated with A 
we have 

A y 0o + [Aq + A %, r) + f(O)]0 o - (0o)r = /i(A o )0o < -Ac0 o , 



or 



A,0 O + [Aq + A (c + 6(j/, r)) + /' (O)]0 O - Mr < 0. (4.4) 
Define u(s,y,r) = e x ° s (f)o(y,T), and we see that u solves 

L c u + f(0)u < 
L c u + f(u) < 0. (4.5) 

To show that c < c*, we argue by contradiction and suppose c* < c. As 
in the construction of the solution (u*,c*), pick 9 > sufficiently small so 
that c* < ce < c < where (ue, cq) solves ()3.3|) . Since eg < c and /e < /, 

L Cg u e + feiue) = 
L c u + < 0. (4.6) 

Replacing u with the function min(-u, 1), we can apply Theorem 2.3 to u = ug 
and u with (3 = eg + b(y, r) and j3 = c + b(y, r) to conclude that for some 
(yo,r ) e T n , f3(y ,To) = P(yo,ro), hence eg = c, a contradiction. Therefore, 
c < c*, which implies 7* < c*. 

Now the opposite inequality 7* > c* will follow from the next lemma, 
analogous to Lemma 6.5 of 2J: 

Lemma 4.1 Let (u,c) be a solution to problem \3. 1)) with u s > 0. Then 
there exists a A > and a function <p(y, r) > such that the function w = 
e Xs (f>(y, t) solves 

L c w + f'(0)w = 0. (4.7) 
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Proof: The proof adapts Lemma 6.5 of Specifically, we define 

0< A = liminf u j s >y> T ) (4.8) 
s^-oo u{s,y,T) 

and we let (sj, yj, Tj) be a sequence such that 

lim " Js >-VrT,) = x (4 g) 

i^oo u{Sj,yj,Tj) 

By the compactness of T n , we may assume (yj, Tj) — > (y^, Too) for some point 
(2/00, Too) G T n . As in the family of functions 

. . m(s + Sj, y,r) , t . 

Wj(s,y,r) = ^- (4.10) 

satisfies Wj > and 

/(u(s + s i; ?/,r) 

L c u> 3 - H — ■ J - —Wj = 0. (4.11) 

u(s + Sj,y,r) 

Parabolic estimates imply that we can take a subsequence of the Wj converg- 
ing in Ci oc (R x T n ) to a function w > that solves 

L c w + f'(0)w = 0. (4.12) 

Also, io(0, |/oo, Too) = 1, so the maximum principle implies w > 0. Now we 
will see that w has the form w = e Xs (p(y, r). Note that for each index j, 

(wj) s (s, y, t) = Us s + s i> y ' T w ( a y r ) (4.13) 
u(a + Sj,y,r) 

so that w s (s,y,r) > Xw(s,y,r) with w s (0, ?/oo, Too) = Aw(0, y^, r^) = A. 
Thus, the function defined by z(s, y,r) = ^ satisfies 2; > A for all (s, y, r) G 
_R x T n with equality holding at the point (0, y^, Too). A straightforward 
computation shows that z solves 



L c z + 2 • V s>y zJ = for (s,y,r) E RxT n . 



(4.14) 



The strong maximum principle then implies that z = A, so that 

d s (we- Xs ) = w s e~ Xs - Xwe~ Xs 

= w (ze- Xs - Xe- Xs ) = 0. (4.15) 
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Hence, the function = w(s,y,r)e~ Xs = 4>(y,r) is independent of s and is 
positive. We conclude that w(s,y,r) = e Xs 4>(y,r). To show that A 7^ 0, we 
plug the function w back into equation ()4.12|) and find that 0(y, r) > solves 

L x cf) = A y (f) + [A 2 + Xb(y, r) + /'(O)]0 - T = -cA0. (4.16) 

From this equation we see that if A = then must be a constant, by the 
maximum principle, hence w must be a constant. However, w cannot be 
constant since ti>(0, 2/00, Too) = 1 and w solves ()4.12|) with /'(0) > 0. This 
proves the lemma. 

Applying this result to the solution (u*,c*), we obtain a function w = 

e Xs 4>(y,r), where 0(y, r) > solves 

L X (p = A,0 + [A 2 + Xb{y, r) + f'(0)]<f> - <f> T = -c*A0, (4.17) 

so that <j)(y, r) is the principal eigenfunction of (J4.1j) with principal eigenvalue 
//(A) = -c*A. Thus 

-c* = ^>inf^ = - 7 *, (4.18) 
A a>o A 

so that c* < 7*. We have now shown that c* =7*. 

Proposition 4.1 For A > ; the function h(X) = attains its infimum 
— c* at a unique point A* > 0. Moreover, h(X) is decreasing in the region 
(0, A*) and is increasing in the region (A*, +00). Thus, no other local minima 
exist. 

Proof: As with (I2.13J) and ()2.14|) . the function h(X) can be expressed as 

W) = ^ = A + W£* (4.19) 
A A J T „ 

where depends on A. The last term is bounded by ||6||oo uniformly in A 
thanks to > 0. From this we see that h(X) — ► 00 as A — > and as A — > +00. 

Now we make use of some results from [7j. Suppose that for some c G R, 
h(X) = —c from some A. Let > be an eigenfunction associated with 
principal eigenvalue /i(A). Then, by (|4.1|) we have 

A y <t> + [A 2 + X(b(y, r) + c) + f (O)]0 - T = (4.20) 
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since c = — Therefore, h(X) = —c if and only if p c (A) = —A 2 — /'(0), 
where p c (A) is the principal eigenvalue and <fi > the principal eigenfunction 
of the operator L defined by 

U = A y( j) + \{b{y, r) + c)<t> -<j> T = p c (\)<p. (4.21) 

According to Lemma 15.2 of [7j, the function A h- > p c (A) is convex. Since the 
function —A 2 — /'(0) is strictly concave, there can be at most two values of 
A such that p c (A) = —A 2 — /'(0). Hence, for any c, there can be at most two 
values of A such that h(X) = —c. This fact, combined with the continuity of 
h(\) and the limits h(X) — > oo as A — > 0, oo, implies the Theorem. 

The variational principle allows the comparison of c* with and without 
an additive time dependent component: 

Proposition 4.2 Writeb(y,r) = b (y)+bi(y, r), whereb (y) = f b(y,r)dr, 
Jq bi(y, t) dr = 0. Let c* sp be the minimal speed when b\ = 0, then \c*\ > \c* \, 
equality holds if b\ = b\ (r) . 

This follows from the analogous comparison property of /i(A), see Theorem 
2.1 of 0. 

5 Variational Principle Based Computation 

The variational formula of Theorem 4.1 is useful from a computational point 
of view because it allows us to compute the minimal speed c* without com- 
puting the solution u of ()1.1|) over long times. The variational formula, on 
the other hand, is a simple minimization problem in just one dimension: for 
a given shear flow, we want to compute 

— c = inf — - — . 

A>0 A 

By Lemma 4.2, the minimum is attained at a unique point A*. For a given 
A > 0, we find p(A) by solving the eigenvalue problem 

A y + [A 2 + Xb(y, r) + /' (O)]0 = /i(A)0. 

We approximate the infimum along the curve A i— > using a Newton's 
method with line search. In this way, our approximation converges quadrat- 
ically in the region near the infimum and decreases with each iteration. Two 
illustrative curves A i— > are shown in Figure two different shears. 
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Because of the periodicity and smoothness of the coefficients, we compute 
the principal eigenvalues //(A) using a Fourier spectral method, as described 
in [2T]. Consider n — 2, y £ [0, 2tt] and r £ [0, 2ir]. For a domain [0, 2n] we 
divide the space into N equally spaced grid points and approximate 

N 

u{y)^^u j S N (y-y j ) (5.1) 

3=1 

for h — ^ and £/j = j/i. The function S^(y) is the band-limited interpolant 
of the Dirac mass 5: 

Sn{v) = ( 2 5i)- (5 - 2) 

Differentiation then translates into matrix multiplication: 

(5.3) 

where is the matrix with entries 

[0 if i = j 

( D *h = | !(-i)C*-i) cot if i ^ J; ( 5 - 4 ) 

and for a second derivative 

= { (5.5) 

Performing a similar discretization in the r dimension, we see that the nu- 
merical eigenvalue problem becomes 

Au = (D%] y + B X - D N>t )u = fi(\)u, (5.6) 

where B\ is the diagonal matrix given by 

5A = (A 2 + f(0))/ + Adiag(6), (5.7) 

and diag(6) denotes the diagonal matrix whose i th diagonal entry is b(y(i), r(i)), 
the value of the function at the point 



y (i) = ((f + 1) mod N)h , r(i) 



% 

N 



h. (5.8) 



In all simulations, we take f(u) = u(l — u). The resulting matrix A is not 
symmetric, since the matrix Dn is not symmetric. We use general double 
precision LAPACK routines [1 to find the princpal eigenvalue, corresponding 
to a positive eigenvector, and thus we obtain points on the curve A i— > ^y-. 
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5.1 Example 1 

First, we define 

b s {y,r) = 5sin(27ry)(1.0 + sin(27rr)), (5.9) 

and we let c*(8) denote the corresponding minimal speed. We consider how 
c*(S) scales with the shear amplitude 5, as in ^7j. It is known that for any 
5 > 0, c*(S) < c* (Recall that in our notation, c* < 0). Supposing c*(5) « 
c*(0) +0(5 P ) we wish to numerically compute the exponent p. According to 
[T7j . we expect c*(S) — Cq is 0(S 2 ) for S small, and 0(5) for 5 large. 

Using the numerically calculated speeds c*(S), we determined the ex- 
ponents p using the least squares method to fit a line to a log-log plot of 
speed versus amplitude. That is, we determined the slope of the best-fit line 
through the data points (log(<5), log(c*(0) — c*(S))) for each shear amplidute 
5. Figure 121 shows this log-log plot. We find that for small amplitude shears 
the minimal speed grows quadratically with the shear amplitude 5: p ~ 1.98. 
For large amplitude shears, the enhancement is linear: p w 1.05, as shown 
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log(8) 

Figure 2: Log-Log plot of enhancement. Slope = 1.98 in the small shear 
amplitude regime. 

next in Example 2. These results agree with but are more accurate and 
computed much faster. 

5.2 Example 2 

Now we consider the effect of temporal frequency on the speed enhancement. 
For each n G {0, 1, 2, 3, 4}, we define 

b s , n (y,r) = 5sin(27ry)(1.0 + sin(27rnr)). (5.10) 

Figure El shows the curves 5 i— »■ —c*(5,n) for each n, plotted simultaneously, 
for small amplitudes. Figure 0] shows the same plot for larger amplitudes. 
Notice that the curves become linear for 5 sufficiently large. When n — 0, the 
shear is independent of time, and the simulations show that the enhancement 
from n = is weakest (cf. Proposition 4.2). For time- dependent shears, the 
enhancment decreases monotonically with increasing temporal frequency. 
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Figure 3: Enhancement by shears with various temporal frequencies, 
quadratic speed growth in the small shear amplitude regime. 

Note that when 5 = 0, bs : n(y,t) = and the minimal speed is c* = —2.0, 
agreeing with the well-known formula 

c* = -2v/f(0) = -2. 
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Figure 4: Enhancement by shears with various temporal frequencies, linear 
speed growth in the large shear amplitude regime. 
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